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Abstract We consider three different models of N non-intersecting Brownian motions on a line 
segment [0,L] with absorbing (model A), periodic (model B) and reflecting (model C) boundary 
conditions. In these three cases we study a properly normalized reunion probability, which, in 
model A, can also be interpreted as the maximal height of N non-intersecting Brownian excursions 
(called "watermelons" with a wall) on the unit time interval. We provide a detailed derivation 
of the exact formula for these reunion probabilities for finite N using a Fermionic path integral 
technique. We then analyse the asymptotic behavior of this reunion probability for large N using 
two complementary techniques: (i) a saddle point analysis of the underlying Coulomb gas and 
(ii) orthogonal polynomial method. These two methods are complementary in the sense that they 
work in two different regimes, respectively for L <C 0(y/W) and L > 0(y/N). A striking feature 
of the large N limit of the reunion probability in the three models is that it exhibits a third- 
order phase transition when the system size L crosses a critical value L = L C (N) ~ \/~N. This 
transition is akin to the Douglas-Kazakov transition in two-dimensional continuum Yang-Mills 
theory. While the central part of the reunion probability, for L ~ L C (N), is described in terms of the 
Tracy- Widom distributions (associated to GOE and GUE depending on the model), the emphasis 
of the present study is on the large deviations of these reunion probabilities, both in the right 
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[L 3> L C (N)] and the left [L <C L C (7V)] tails. In particular, for model B, we find that the matching 
between the different regimes corresponding to typical L ~ L C (N) and atypical fluctuations in 
the right tail L 3> L C (N) is rather unconventional, compared to the usual behavior found for the 
distribution of the largest eigenvalue of GUE random matrices. This paper is an extended version of 
[G. Schehr, S. N. Majumdar, A. Comtet, J. Randon-Furling, Phys. Rev. Lett. 101, 150601 (2008)] 
and [P. J. Forrester, S. N. Majumdar, G. Schehr, Nucl. Phys. B 844, 500-526 (2011)]. 

1 Introduction 

Non-intersecting random walkers, first introduced by de Gennes [T], followed by Fisher [2] (who 
called them "vicious walkers"), have been studied extensively in statistical physics as they appear 
in a variety of physical contexts ranging from wetting and melting all the way to polymers and 
vortex lines in superconductors. Lattice versions of such walkers also have beautiful combinatorial 
properties [3j. Non-intersecting Brownian motions, defined in continuous space and time, have also 
recently appeared in a number of contexts. In particular their connection to the random matrix 
theory has been noted in a variety of situations [4.5.6.7.8.9,10.1 lirT2 | ri3| J T4 ,15.16. 17ll 18| . Rather 
recently, we have unveiled an unexpected connection between vicious walkers and two-dimensional 
continuum Yang-Mills theory on the sphere with a given gauge group G [19) . which depends on 
the boundary conditions in the vicious walkers problem (a different type of connection between the 
Yang-Mills theory and vicious walkers problems had also been noticed in Ref. |20jV 

Specifically, we consider a set of N non-intersecting Brownian motions on a finite segment [0, L\ 
of the real line with different boundary conditions. Assuming that all the walkers start from the 
vicinity of the origin, we then define the reunion probability as the probability that the walkers 
reunite at the origin after a fixed interval of time which can be set to unity without any loss of 
generality. Within the time interval [0,1], the walkers stay non- intersecting. Next we 'normalize' 
this reunion probability in a precise way to be defined shortly. In one case, namely when both 
boundaries at and L are absorbing, one can relate this 'normalized' reunion probability to the 
probability distribution of the maximal height of N non-intersecting Brownian excursions. In Ref. 
|19| it was shown that this normalized reunion probability in the Brownian motion models maps 
onto the exactly solvable partition function (up to a multiplicative factor) of two-dimensional Yang- 
Mills theory on a sphere. The boundary conditions at the edges and L select the gauge group 
G of the associated Yang-Mills theory. We consider three different boundary conditions: absorbing 
(model A), periodic (model B), and reflecting (model C) which correspond respectively to the 
following gauge groups in the Yang-Mills theory: (A) absorbing — s- Sp(27V), (B) periodic — > U(iV) 
and (C) reflecting — > SO(2iV). As a consequence of this connection, in each of these Brownian 
motion models, as one varies the system size L, a third order phase transition occurs at a critical 
value L = L C (N) = 0{\J~N) in the large N limit. It was shown in Ref. [21] that a similar third 
order phase transition also occurs for the probability distribution function of the largest eigenvalue 
of random matrices belonging to the standard Gaussian ensembles. Furthermore, third order phase 
transitions in the large deviation function of appropriate variables have also been found recently in 
a variety of other problems, such as in the distribution of conductance through a mesoscopic cavity 
such as a quantum dot [22,23,24 and in the distribution of the entanglement entropy in a bipartite 
random pure state |25U26) . 

Close to the critical point, these reunion probabilities in the Brownian motion models, properly 
shifted and scaled, can be related to the Tracy- Widom (TW) distributions. Let us briefly remind 
the readers about the Tracy- Widom distributions. Tracy- Widom distribution describes the limiting 
form of the distribution of the scaled largest eigenvalue in the three classical Gaussian random 
matrix ensembles |27| . These limiting distributions are usually denoted by Fi{t) for the Gaussian 
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orthogonal ensemble (GOE), by ^(t) for the Gaussian unitary ensemble (GUE) and by Fn{t) for 
the Gaussian symplectic ensemble (GSE). For example, in the GUE case where one considers the 
set of TV x TV complex Hermitian matrices X with measure proportional to e~ TlX and denote A max 
the largest eigenvalue, the typical fluctuations around its mean value (A max ) — V2N have a limiting 
distribution [2"T] 

Jin^Pr^V^A^ - V2N) < tj = exp J (s - t)q 2 (s)dsj := F 2 (t), (1) 

known as the /3 = 2 Tracy- Widom distribution. In Eq. (TTJ), g(s) satisfies Painleve II (PII) differential 
equation 



with the asymptotic behavior 



q"( S ) = S q(s) + 2q 3 ( S ) , (2) 



q(s) ~ Ai(s) as s — ?> oo , (3) 



where Ai(s) is the Airy function. One can indeed show that this asymptotic behavior ([3]) determines 
a unique solution of PII @, known as the Hastings-McLeod solution. Similarly, the distribution 
T\ concerns real symmetric matrices. Explicitly, with the GOE specified as the set of TV x TV 
real symmetric matrices X with measure proportional to e~ TrX / 2 , and A max denoting the largest 
eigenvalue, one has [28] 

lirn^ Pr ^ V2TV5 (A max - V2N) < t J = exp ( - i jf ((« - t) g 2 (s) + ?(«)) 

:=Ji(t), (4) 

known as the /3 = 1 Tracy- Widom distribution. Because of their relevance in many fundamental 
problems in mathematics and physics, these TW distributions have been widely studied in the 
literature. In particular, their asymptotic behaviors are given, to leading order by 

>/j(t) ~exp(-^|T| 3 ) , i->-oo, 

(5) 

I - Tp(t) - exp (-^t 3/2 ) , M oo 

where /? = 1 and (3 = 2 correspond respectively to the GOE and the GUE case. 

Our study of constrained vicious walkers problems started in Ref. |12) where we derived, using a 
Fermionic path integral method, an exact expression for the ratio of reunion probability, for model 
A and for any finite number TV of walkers. This calculation, in this specific case, was motivated by 
the interpretation of this ratio in terms of an extreme value quantity, namely the maximal height of 
TV non- intersecting Brownian excursions, which is recalled below in section [2] (see also Fig. [I). We 
showed later in Ref. |19| that this ratio, and its extension to other boundary conditions in model 
B and C, is actually equal, up to a multiplicative prefactor, to the partition function of Yang-Mills 
theory on the sphere with a given gauge group G, which depends on the boundary conditions in 
the vicious walkers problem (see Table [T]). Following the pioneering works of Refs. |291l30] where 
the large TV analysis of Yang-Mills theory on the sphere with the gauge group G = U(TV), and 
extended to other classical Lie groups, including G = Sp(2TV) or G = SO(2TV), in Ref. [31], we 
could also perform in Ref. [19] the large TV analysis of these reunion probabilities. However, most of 
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these results were announced in Refs. [T2"lll9j without any detail. The purpose of the present paper 
is to give a self-contained derivation of these results, both for finite and large N. For large N we 
characterize not only the distribution of typical fluctuations, which can be expressed in terms of T\ 
and Ti but also provide a detailed analysis of the large deviations of these reunion probabilities, 
characterizing atypical fluctuations. A special emphasis is put on the matching between different 
regimes (typical and atypical) of the reunion probability as a function of L. In particular, for 
the case of periodic boundary conditions, we find that the matching found in this vicious walker 
problem is quite different from the corresponding matching observed in the distribution of the 
largest eigenvalue of Gaussian random matrices in the unitary ensemble |I 31f32 , 33 , 34 . 35 , 551137) . 

The paper is organized as follows. In section 2 we describe the three models A, B and C and 
summarize the main results for the ratio of reunion probabilities in the large N limit. In section 3 
we give the details of the path integral method which allows us to obtain an exact expression of 
these ratios in each of the three models for any finite N and L. In section 4, we compute the large 
N limit of this ratio for model A and analyze in detail the typical fluctuations (namely the central 
part of the distribution) as well as the large deviations: both for the left and for the right tails. 
The corresponding large N analysis for model B and C are performed in section 5 and section 6, 
respectively, before we conclude in section 7. Some details have been relegated to Appendix A. 



2 Models and main results 

We consider three different models of N non-intersecting Brownian walkers on a one-dimensional 
line segment [0, L] and label their positions at time r by x\(t) < . . . < xn(t). These three models, 
denoted by model A, B and C, differ by the boundary conditions which are imposed at x = and 
x = L: 

• in model A, we consider absorbing boundary conditions both at x = and x = L, 

• in model B, we study periodic boundary conditions, which amounts to consider N non-intersecting 
Brownian motions on a circle of radius L/2ir, 

• and in model C we consider reflecting boundary conditions both at x — and x = L. 





Model A 


Model B 


Model C 


Boundary conditions at x = and x = L 


absorbing 


periodic 


reflecting 


Corresponding affine Weyl chamber 


Cat 


A N -! 


Bn 


Gauge group of the associated YM2 theory 


Sp(2iV) 


X1(N) 


SO(2iV) 



Table 1 Summary of our results for the three different models A, B and C. 



Model A: In the first model the domain is the line segment [0, L] with absorbing boundary condi- 
tions at both boundaries and L. This corresponds to A-dimensional Brownian motion in an affine 
Weyl chamber of type CV |39U40| . The N non- intersecting Brownian motions start initially at the 
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positions, say, {ei, ei, ■ ■ ■ , €n} m the vicinity of the origin where eventually we will take the limit 
€i — > for all i, as shown later. We define the reunion probability i?^(l), where the superscript 'A' 
refers to model A, as the probability that the walkers return to their initial positions after a fixed 
time t = 1 (staying non-intersecting over the time interval r <E [0,1]). We define the normalized 
reunion probability 

ML) = fg, (6) 

such that liniL-j.00 Fm{L) = 1. This ratio becomes independent of the starting positions e^'s in the 
limit when e; — > for all i, as shown later. Hence, Fn(L) depends only on N and L. 

This ratio Fm{L) in Model A has also a different probabilistic interpretation. Consider the same 
model but now on the semi-infinite line [0,oo) with still absorbing boundary condition at 0. The 
walkers, as usual, start in the vicinity of the origin and are conditioned to return to the origin exactly 
at r = 1 (see Fig. |T|) . If one plots the space-time trajectories of the walkers, a typical configuration 
looks like half of a watermelon (see Fig. [TJ, or a watermelon in presence of an absorbing wall. 
Such configurations of Brownian motions are known as non-intersecting Brownian excursions |10] 
and their statistical properties have been studied quite extensively in the recent past. A particular 
observable that has generated some recent interests is the so-called 'height' of the watermelon \12\ 
IM l l4T | l42 i r43 i r44 l l45 p 46| defined as follows (see also Ref. [47] for a related quantity in the context of 
Dyson's Brownian motion). Let Hn denote the maximal displacement of the rightmost walker xn 
in this time interval r € [0,1], i.e., the maximal height of the topmost path in half-watermelon 
configuration (see Fig. Q]), i.e., H N = max r {xjv(i"), < r < 1}. This global maximal height H N 
is a random variable which fluctuates from one configuration of half-watermelon to another. What 
is the probability distribution of iijv? For N = 1 the distribution of Hpj is easy to compute and 
already for N = 2 it is somewhat complicated [42 . In Ref. [T^] an exact formula for the distribution 
of iJ/v, valid for all N, was derived using Fermionic path integral method, which we remind below 
for consistency (see also [42" l l43 [ l44] for a derivation using different methods). The distribution of 
-ffjVi m the large N limit, is quite interesting as it gives, in the proper scaling limit, the distribution 
of the maximum (on the real line) of the Airy2 process minus a parabola [6,48,49 . This latter 
process describes the universality class of the Kardar-Parisi-Zhang (KPZ) equation in the so-called 
"droplet" geometry. It was known rather indirectly from the work of Ref. [6] on the Airy2 process 
that the limiting distribution of Hn should then be given by T\ . This was one of the main result of 
Ref. |19| to show this result by a direct computation of the limiting distribution of i/jv- This result, 
for the vicious walkers problem, was then recently proved rigorously in |46| using Riemann-Hilbert 
techniques. In Ref. [50] a direct and rigorous proof was established that the distribution of the 
maximum of Airy2 minus a parabola is indeed given by T\. We refer the reader to Refs. [51,52,53] 
for more recent results on the extreme statistics of the Airy2 process minus a parabola. We also 
mention that the limiting distribution of Hm was measured in recent experiments on liquid crystals, 
and a very good agreement with J-\ was indeed found [54] . 

To relate the distribution of Hn in the semi-infinite system defined above to the ratio of reunion 
probabilities in the finite segment [0, L] defined in ((6]), it is useful to consider the cumulative 
probability Pr(iJ/v < L) in the semi-infinite geometry, where L now is just a variable. To compute 
this cumulative probability, we need to calculate the fraction of half- watermelon configurations (out 
of all possible half- watermelon configurations) that never cross the level L, i.e., whose heights stay 
below L over the time interval r e [0, 1] (see Fig. Q]). This fraction can be computed by putting an 
absorbing boundary at L (thus killing all configurations that touch/cross the level L). It is then 
clear that Pt(Hn < L) is nothing but the normalized reunion probability Fn(L) defined in ([5]). 
As mentioned above, this cumulative probability distribution of the maximal height was computed 
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exactly in Ref. JT2] 

F N (L) :=Pt(H n <L) 

.00 00 
An \ ~* 



£2N 2 +N ' ' ' 

ni — — oc un — — c 



IV 2 

Zi 2 (n?,...,n^)(nn|) e -^^^, (7) 



where 



Av(i/i,...,2/jv) = {Vi-Vj) (8) 

l<i<j<N 

is the Vandermonde determinant and where the amplitude is given by 

n 2N 2 +N 

An = m — "1 ■ (9) 

2 A-+A72 U N-1 r(2 + i)r(3/2 + ?) 

This quantity Fn{L) was also computed in Refs. [4U1FH] using different methods. In Ref. |43j 
the authors studied the maximal height of N non-intersecting Brownian excursions and used the 
Karlin-McGregor formula [55]. In Ref. [44] the author computed the cumulative distribution of the 
maximal height of A^-non intersecting discrete lattice paths (of n discrete steps) in presence of a 
wall, using the Lindstrom-Gessel-Viennot formula |56U57) . and then he considered the asymptotic 
limit of a large number of steps n. The formulas obtained using these two similar methods are 
actually different from the one given in Eq. ([7]) obtained by using a fermionic path integral and it 
is a non-trivial task to check, as expected, that they are indeed equivalent to ([7]) [43] . A derivation 
of the formula given above in Eqs. ([71 ©j which turns out to be the most convenient expression of 
Fn{L) in view of a large N asymptotic analysis, is given in section[3l Remarkably, if one denotes by 

Z(A,G) (10) 

the partition function of the two-dimensional (continuum) Yang-Mills theory on the sphere (denoted 
as YAI2) with gauge group G and area A it was shown in Ref. [19] that Fn{L) is related to YM2 
with the gauge group G = Sp(2A^) via the relation 

F n (L)ocz(a=^-N,S P (2N)^ . (11) 

In Ref. [29113 1], it was shown that for large N, Z(A, Sp(2A^)) exhibits a third order phase transition at 
the critical value A = ir 2 separating a weak coupling regime for A < ir 2 and a strong coupling regime 
for A > tt 2 . This is the so called Douglas-Kazakov phase transition [29], which is the counterpart 
in continuum space-time, of the Gross- Witten-Wadia transition [58,59J found in two-dimensional 
lattice quantum chromo-dynamics (QCD), which is also of third order. Using the correspondence 
L 2 = 2ir 2 N/A, we then find that Fn(L), considered as a function of L with N large but fixed, also 
exhibits a third order phase transition at the critical value L C (N) = s/2N. Furthermore, the weak 
coupling regime (A < tt 2 ) corresponds to L > \/2N and thus describes the right tail of Fn(L), 
while the strong coupling regime corresponds to L < \/2N and describes instead the left tail of 
Fn(L) (see Fig. The critical regime around A = tt 2 is the so called "double scaling" limit in 
the matrix model and has width of order N~ 2 / 3 . It corresponds to the region of width 0(N^ 1 ^ 6 ) 
around L = \/2N where Fn(L) : correctly shifted and scaled, is described by the Tracy- Widom 
distribution J"i(t) in Eq. ([4]). 
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Fig. 1 Trajectories of N non-intersecting Brownian motions x\(t) < 2:2(1") < . . . < xn(t), all start at the 
origin and return to the origin at r = 1, staying positive in between. Fn(L) denotes the probability that the 
maximal height Hn = max T {ijv(r),0 < r < 1} stays below the level L over the time interval < t < 1. 



Although the occurrence of the Painleve transcendent q(s) © in this double scaling limit was 
known since the work of Periwal and Shevitz |60| . its probabilistic interpretation in relation to the 
Tracy- Widom distribution was one of the main achievements of Ref. |19] . In this paper, we provide 
a detailed analysis of the three regimes: the left tail, the central part and the right tail of Fn(L). 
Our results can be summarized as follows 

'F N (L) ~ exp [-N 2 ^ (L/VW)] , L < V2N & \L- V2N\ ~ 0(VN) 
< F N (L) - T x [2 u /6at1/6( L _ 727V)] j L ~ y/2N k \L - V2N\ - ©(A^- 1 / 6 ) (12) 

^ 1 - F N (L) - exp [-N<fr$ (L/Vw)] , L > V2N & \L - V2A| - 0(VN) , 

where J-\ is the TW distribution for GOE, whose explicit expression is given in ([4]) and its asymptotic 
behaviors are given in Eq. ((SJ). The rate functions 4>±{x) can be computed exactly (see later): of 
particular interest are their asymptotic behaviors when L — > \/2N from below (left tail) and from 
above (right tail), which are given by 

1 fi 

^(X)~y(l-X) 3 ,x^r, (13) 
o9/2 
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Fig. 2 Left: Schematic sketch of Fn(L) as defined in Eq. (0 for N vicious walkers on the line segment [0, L] 
with absorbing boundary conditions at both ends, as a function of L, for fixed but large N. Right: Sketch 
of the phase diagram in the plane (A, N) of two-dimensional Yang-Mills theory on a sphere with the gauge 
group Sp(2iV) as obtained in Ref. [291131) . The weak (strong) coupling phase in the right panel corresponds 

to the right (left) tail of F^(L) in the left panel. The critical region around A — tv 2 in the right panel 

corresponds to the Tracy- Widom (TW) regime in the left panel around the critical point L C (N) = \/2N. 



The different behavior of F N {L) in Eq. {T^J for L < V2N and L > V2N leads, in the limit N -> oo 
to a pha se transition at the critical poi nt L = \/2N in the following sense. Indeed if one scales L 
by y/2N, keeping the ratio x = Lj\j2N fixed, and take the limit N — > oo one obtains 



lim 

N—>oc 



1 



■]nF: 



N X 



<t>-{x) , X < 1 
, x > 1 . 



(14) 



If one interprets Fm{L) in Eq. ([7]) as the partition function of a discrete Coulomb gas, its logarithm 
can be interpreted as its free energy. Since 4>^_(x) ~ (1 — x) 3 when x approaches 1 from below, then 
the third derivative of the free energy at the critical point x = 1 is discontinuous, which can then 
be interpreted as a third order phase transition. 

On the other hand, comparing the asymptotic behavior of T\ in Eq. ([5]) with the ones of the rate 
functions (|13j) we can check that the expansion of the large deviation functions around the transition 
point coincides with the tail behaviors of the central region Tracy- Widom scaling function. This 
property holds here both fo r th e left and the ri ght tails. For instance, consider first the left tail in 
Eq. (|12p. i.e., when L <C \/2N. When L — > \2N from below, we can substitute the asymptotic 
behavior of the rate function (x) from Eq. (fl3|) in the first line of Eq. (fT2| . This gives 



F N (L) - exp 



2 5/2 



N 1/2 {V2N - L) 3 ) , l<V27V-£< 



(15) 



On the other hand, consider now the second line of Eq. (TT2|) that describes the central typical 
fluctuations. When the deviation from the mean is large, i.e., s/2N— L ~ 0{\f~N) (compared to the 
typical scale - 0(N~ 1 / & )). we can substitute in the second line of Eq. ([HZ]) the left tail asymptotic 
behavior of the Tracy- Widom function J-\(t) as described in the first line of Eq. ((5]) (with (3 = 1). 
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This gives, 



F N (L) ~ exp 



~{2 n l 6 N x '\V2N-L)J 



(16) 



which, after a trivial rearrangement, is identical to the expression in Eq. (|15[) . Thus the left tail of 
the central region matches smoothly with the left large deviation function. Similarly, on the right 
side, using the behavior of 4>+(x) in Eq. (|13p. one finds from Eq. (jT2"j) . that 

1 - F N (L) ~ exp f — N 1/4 {L - V2N) 3/2 J , 1«L - V2N < V2N , (17) 

which matches perfectly with the right tail of the central part described by JF\ (i) (0 [T2| . Hence 
in this case of model A, the matching between the different regimes is similar to the one found 
in previous studies of large deviation formulas associated with the largest eigenvalue of random 
matrices [T51I521I551IM1I551I5S] . 

Model B: In the second model we consider periodic boundary conditions on the line segment [0, L\. 
Alternatively, one can think of the domain as a circle of circumference L (of radius L/2n). This 
corresponds to AT-dimensional Brownian motion in an affine Weyl chamber of type A n _i [39 , TO] . All 
walkers start initially in the vicinity of a point on the circle which we call the origin. We can label 
the positions of the walkers by their angles {#1, 82, ■ ■ ■ , On} (see section [3] for details). Let the initial 
angles be denoted by {ei, £2, ■ • ■ , ejv} where e^'s are small. Eventually we will take the limit ei — > 0. 
We denote by Rf(l) the reunion probability after time r — 1 (note that the walkers, in a bunch, 
may wind the circle multiple times), i.e, the probability that the walkers return to their initial 
positions after time r = 1 (staying non-intersecting over the time interval r 6 [0,1]). Evidently 
depends on N and also on the starting angles {ei, €2, ■ ■ ■ , £iv}- To avoid this additional 
dependence on the e^'s, let us introduce the normalized reunion probability defined as the ratio 

om = jgw , (is) 

where we assume that we have taken the ej — > limit. One can actually check that the ej's depen- 
dence actually cancels out between the numerator and the denominator: this is the main motivation 
for studying this ratio of reunion probabilities (fT8)) in this case. Although in the case of model A it 
is rather natural to expect that this limit q — > is well defined, given the interpretation of Fjy(L) 
as the cumulative distribution of the maximal height of non-intersecting excursions, this is not so 
obvious for Gm{L) where such an interpretation does not exist. Nevertheless one can show that this 
property also holds in this case, yielding (see Ref. Q25] and also section [3]) 

Gn(L) = ^ £ •■• E A%(m,...,n N )e- 2 -£^^ , (19) 



where An(jii, ■ ■ ■ ,tin) is the Vandermonde determinant ([8J and the prefactor 
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ensures that Gn{L — > oo) = 1. In Ref. [19] it was shown that this normalized reunion probability 
Gn{L) is, up to a prefactor, exactly identical to the partition function of the 2-d Yang-Mills theory 
on a sphere with gauge group \J(N): 



4ir 2 N \ 
A=— ,U(JV) 



(21) 



This partition function Z (A,XJ(N)) exhibits the Douglas-Kazakov third order phase transition for 
A = 7r 2 which means that in that case the transition between the right tail and the left tail behavior 
of Gn(L) occurs for L = L C (N) = VAN. In this paper, we provide a detailed analysis of the various 
regimes, right tail, central part and left tail of Gn(L). Our results can be summarized as follows 



f G N (L) 



exp [~N 2 <f)B (l/Van)] , L < Van k\L- 



G N {L) ~ J" 2 [2 2 / 3 A fl / 6 (i — \/4JV)] , L ~ VAN k.\L — V~A~N\ ~ 0(N~ 1 / 6 ) (22) 

1 - G N {L) ~ (-1) N exp [-JV0f (L/VAN)] , L > VAN & \L - V^N\ - 0{V~N) , 

where Ti is the TW distribution for GUE, whose explicit expression is given in (JlJ and its asymptotic 
behaviors are given in Eq. ([5]). The rate functions 4>±(x) = \4>±{x) can be computed exactly (see 



later). Of particular interest are their asymptotic behaviors when L 
and from above (right tail), which are given by 



MJV from below (left tail) 



B (?) 



-(1 - x 3 , X 
3 V ' ' 

2 7/2 



1 



(23) 



3/2 



, X 



Notice also the oscillating sign in the right tail of Gn{L) in Eq. (j2"2")l which is not problematic here 
as Gm{L) does not have the meaning of a cumulative distribution. As explained before for model 
A in Eq. (|14p . the cubic behavior of <j>^.(x) when x approaches 1 from below is again the signature 
of a third order phase transition in this model. On the other hand, by comparing the asymptotic 
behavior of T<z in Eq. ([S]) with the one of the rate functions (|2"3"|) we can check, in this case of 
model B, that the expansion of the large deviation functions around the transition point L = V AN 
coincides with the tail behavior of the TW scaling function, only for the left tail. 
Indeed, using the behavior of 4>-{x) in Eq. (|23|) . one finds from Eq. ([22]) that 



Gn(L) ~ exp 



1n 1 /\l- Van) 3 



i < Van - l < 



(24) 



which matches perfectly with the left tail of the central region ([5]). However, we find that this 
property does not hold for the right tail. Indeed, using the asymptotic behavior of </>5 (x) in Eq. (|23|) 
one finds from (l22l) 



1-Gjv(L)~ (-irexp 



N^(L 



1 < L 



(25) 
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double scaling 


large deviation 




regime 


regime 




e>(7v-i/6) 


L > V^N 




1 ► 

-Across L 



Fig. 3 Crossover from the double scaling to the large deviation regime which occurs only for model B. 
This crossover happens for L cross ~ iV~ 1/ ' 6 (lniV) 2/ ' 3 . 



with an oscillating sign oc (— 1) . On the other hand the right tail of the central region, described 
by F 2 {t) in Eq. ©, yields for L - y/AN > A -1 / 6 

l-G N {L)~l-F 2 

without any oscillating sign and where the argument of the exponential is twice larger than the 
one in (|25[) . This mismatching is the sign of an interesting crossover which happens in this case 
and which we study in detail below. It can be summarized as follows. Close to L = s/4N, with 
L — y/AN ~ 0(N~ x / & ), a careful computation beyond leading order shows that 

lnGjv(L) = In F 2 (t) + (-l) JV - 1 JV- 1 / 3 2 _1 / 3 ?(t) + 0(N~ 2 / 3 ) , t = 2 2 / 3 A 1 / 6 (L - y/AN) , (27) 

where q(t) is the Hastings-McLeod solution of PII In the large t limit, these two competing 

terms In F 2 {t) and q(t) in Eq. behave like ([31 EJ) 

| In Ja(*)| ~ e--^' 2 , q(t) ~ Ai(t) ~ . (28) 

Therefore what happens when one increases L from the critical region L-ViN ~ 0(N~ 1 / & ) towards 
the large deviation regime in the right tail L > V 'AN is the following (see Fig. [3]): the amplitude of 
the second term in the right hand side of Eq. (|2T|) which is oscillating with N, increases relatively 
to the amplitude of the first term. At some crossover value L = L cross (N) it becomes larger than 
the first one and in the large deviation regime it becomes the leading term, still oscillating with N 
(22) . Balancing these two terms and making use of the above asymptotic behaviors (|28p one obtains 
an estimate of L cmss (N) as 

L cmas - VAN ~ N-^iln A0 2/3 . (29) 

Note that such a peculiar crossover is absent in the distribution of the largest eigenvalue of GUE 
random matrices and it is thus a specific feature of this vicious walkers problem. 

Model C: We consider a third model of non- intersecting Brownian motions where the walkers 
move again on a finite line segment [0,L], but this time with reflecting boundary conditions at 
both boundaries and L. This corresponds to ^-dimensional Brownian motion in an affine Weyl 
chamber of type Bjy |39lH0] . Again the walkers start in the vicinity of the origin at time r = and 



2 2/3 N l/6( L , 



exp 



(26) 
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we consider the reunion probability i?£(l) that they reunite at time r = 1 at the origin. Following 
Models A and B, we define the normalized reunion probability 



E N {L) 



(30) 



that is independent of the starting positions {ex, 62, ... , ejv} in the limit when all the e-s's tend to 
zero and hence depends only on N and L. As shown in Ref. [IS], see also in section [31 En(L) can 
be computed exactly as 



E N {L) = 



L 2N 



2 -N 



2-ij = \ 



(31) 



ni— — oo riN — — oo 



where An(i)i, ■ ■ ■ , j/jv) is the Vandermonde determinant (|58[) and the prefactor Cjv is given by 

2N 2 -N 2N/2-N 2 



c 



N 



nf=o lr (2+i)r(i/2 + i) 



(32) 



which ensures that Ejy(L — > 00) = 1. In Ref. PH] we showed that En(L), up to a prefactor, is 
exactly identical to the partition function of the 2-d Yang-Mills theory on a sphere with gauge 
group SO (2 AO: 



E N (L) oc Z f A = ^- N, SO(2A) J 



(33) 



This partition function Z (A,SO(2N)) exhibits the Douglas-Kazakov third order phase transition 
for A = tt 2 which means that in that case t he tr ansition between the right tail and the left tail 
behavior of En(L) occurs for L = L C {N) = y/2N. In this paper, we provide a detailed analysis of 
the various regimes, right tail, central part and left tail of En(L). Our results can be summarized 
as follows 



(E N (L) 



E N (L) 



[l-E N (L) 



exp [-N 2 <j) c (L/V2N)] , L < V2N k\L- ^/2A| - 0{yfN) 



Ms) 



s = 2 11 / 6 A 1 / 6 (L - V2N) , L 



/2A & \L - y/2N\ - 0{N- 1 / Q ) (34) 



exp 



[-N<j><l (L/V2N)] , L > V2N k\L- V2N\ - O(VN) , 



where T2 is the TW distribution for GUE and T\ the TW distribution for GOE, whose explicit 
expressions are given in (JT1 |4j) and their asymptotic behaviors are given in Eq. ((5]). The rate functions 
4>±(x) = <fi±(x) can be computed exactly (see below): their asymptotic behaviors are given in Eq. 
(fT3| . As explained before for model A in Eq. (fT4|). the cubic behavior of 4>_(x) when x approaches 
1 from below is again the signature of a third order phase transition in this model. Similarly, as in 
the case of model A, we can check that the expansion of the large deviation functions around the 
transition point L = \/2N coincides with the tail behaviors of the central region both in the left 
tail 



E N {L) 



■ exp 



2 5/2 



N 1/2 (V2N- Lfj 



1 < V2A - L < V2N 



(35) 
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as well as in the right tail 

/ 2 15 / 4 \ 

1 - E N (L) ~ - cxp ( — TV 1/4 (L - V2N) 3/2 J , 1< L - V2N < V2TV , (36) 

where, here again, the minus sign is not problematic as Ej^(L) does not have the interpretation of 
a cumulative probability distribution. 



3 Derivation of the formula for the reunion probabilities 

In this section, we derive the expressions of the normalized reunion probabilities given in Eqs. ([JJ 
ri"T)l l3"Tj) . The derivations are based on a Fermionic path integral method. For model A, this result 
was first reported in Ref. [12] . 



3.1 Model A: absorbing boundary conditions at x = and x = L 
We start by the computation of F^{L) defined in Eq. ©. It then follows that 

N{e,L) 



F N {L) = lim 



TV(e, L — » 00) 



(37) 



where we use the notation e = ei, • • • , ejv and where TV(e, L) is the probability that the TV Brownian 
paths, with diffusion coefficient D = 1/2, starting at < c\ < ■ ■ ■ < ejv at r = come back to 
the same points at r = 1 without crossing each other and staying within the interval [0,L], with 
absorbing boundary conditions at x = and x = L. 

To proceed, let us first consider the simple case of TV independent and free Brownian walkers 
on a line, each with diffusion constant D = 1/2, over the unit time interval [0, 1], but without the 
non-intersection constraint. The probability measure of an assembly of TV trajectories, that start 
at e and also end at e, would then be simply proportional to the propagator 



Pf ree (e, T = 0|e,T = 1) 



T>Xi (r) exp 



1 N 

i=l 



dxi 
d7 



dr 



N 



~\5[x i (Q)-e i ]8[x i (l)-e i ] . (38) 



If, in addition, we subject the walkers to stay within the box [0,L] during the time interval [0, 1], 
this is equivalent to putting an infinite potential at the two ends of the box [0, L\. Then one can 
use path integral techniques to write TV(e, L) as the propagator 



N(e,L) = (e\e-" L \e) 
where Vl{x) is a confining potential with 

V L (x) = 



1 d 2 
-2dx? + VL{Xi) 



(39) 



, x G [0, L] , 
+oo , x <£ [0, L] 



(40) 
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If one denotes by E the eigenvalues of Hl ([59)1 and \E) the corresponding eigenvectors, one has 



N(e, L) = \ |^(e)| 2 e-^ , V E (e) = (e\E) 



E 



(41) 



So far, we have not implemented the non-intersection constraint. The important observation 
that we make is that this constraint can be incorporated within the path integral framework by 
simply insisting that the many body wave function &e(x) = We{x\, ■ ■ ■ , xjv) must be Fermionic, 
i.e. it vanishes if any of the two coordinates are equal. This many-body antisymmetric wave func- 
tion is thus constructed from the one-body eigenfunctions of Hl by forming the associated Slater 
determinant. In the case of model A, we note that the single particle wave-functions vanishing at 
x = and x = L are 



<f> n (x) = J ^ sin {^jpj , neZ + , 

such that the eigenfunctions ^(e) and eigenvalues E in Eq. (j4"Tjl are given here by 

1 



<Me) = 



det <f>m(ti) i = 777^7 y^"-? 

/]VT l<i,j<N Y ' y 31 2L 2 ^ 1 

2—1 



(42) 



(43) 



Hence one has 



N(e,L) 



1 

N\ V L 



N OO 



E-E 

m— l njv — l 



sin 


' ni7rei ^ 
> L < 


sin 


' ni7re 2 ^ 


■ ■ ■ sin 


' ni-irejv \ 


sin 


' n 2 7rei ' 
v L / 


sin 


v L i 


• • • sin 


> n 2 TTe N j 



njy7rei 
v L , 



n N ire 2 
> L , 



' ri W 7rejy 
v L 



e ai- 



(44) 



We now have to study the limit of N(e, L) in Eq. (|4"4")l when e 2 ; — >• (Vi = 1, • • • , N). One can then 
check that to leading order, one has 



sin 


' ni7rei > 
> L < 


sin 


' Til 7T£2 ^ 
> L < 


■ ■ ■ sin 




sin 


' n 2 7rei > 


sin 


' n 2 7re 2 < 
k L / 


• • ■ sin 


> n 2 ire JV j 



sin(2^pL) sin(^ 



n N Tte N 
L , 



(l L 



2JV 2 



2\2 



AT 



e 2 - e 2 ) 2 



<<j 



where a\ is a numerical constant independent of n^s and e^'s. Therefore one obtains 



(45) 



F N {L) 



"2 



iV(e, L — » oo) 



£ 2 ) 2 



JV 



Ef =1 < 



n i — 1 



riiv — 1 i— 1 



2<J 



(46) 

with a2 = 2 Oi. To compute N(e, L — > oo), starting from Eq. (|4~4"|) one first notices that the product 
of determinants can be replaced by the same limiting behavior as in the e$ — > limit given in Eq. 
(|45| as this product is a function of the variables e^/L's only. One can then perform the remaining 
multiple sums over n^s in f|44[) by noticing that they can be replaced by integrals in the limit 
L —> oo - as done below in Eq. (f5U|) . Finally, one obtains that 



lim 



nLe?rw^-ef) 2 

N(e, L — > oo) 



= a 3 , 



(47) 
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where 03 is a constant independent of L. Therefore, combining these results ([4l)I |4"T)) , and using the 
symmetry of the summand in Eq. (|46|) under the transformation rii — > — rij, one arrives at 



, 00 00 N 

w)=j0tw e ••• e n^ 2 -^-^"* m 

ni— — 00 n^ — — 00 i— 1 i<ij 

as announced before (O, where the constant An remains however undetermined. It can be computed 
using the normalization condition of Fn(L) which follows directly from its definition ^ 

lim F N (L) = 1 . (49) 

L— >oo 

Indeed, as we mentioned it above, in the limit when L — > 00, the discrete variables im.i/ L = ki which 
naturally enter into the expression of Fn(L) in Eq. (|48j) become continuous variables. Therefore 
the discrete sums in Eq. (|48[) can be replaced by integrals in the L — > 00 limit, giving 



L 



lim F N (L) = / dh--. dk N ]J kf - tf?^ ^ k > . (50) 

n J ~°° i=i i<j 



If one restricts the integrals over ki E [0, +00) and performs the change of variable Xi = kf /2 one 
obtains 



(51) 



^ 1 i<j 

where the integral can now be evaluated using a limiting case of the Selberg integral [64J: 

•fax... f ^ n - n f ^ = n r(1 + Vn+T + - ■ < 52 > 

i<j i=l j=0 *> ^ 



JO 



yielding (with 7=1 and a = 3/2). the formula for An given in Eq. [see also Eq. 



3.2 Model B: periodic boundary conditions 

Here we consider N non- intersecting Brownian motions with diffusion coefficient D = 1/2 on a line 
segment [0, L] with periodic boundary conditions or equivalently on a circle of radius L/2n, starting 
in the vicinity of the origin at time t = 0. In this case, the computation of the ratio of reunion 
probability in Eq. (fT9| can be done along the same line as before, (|39ti43|) . In this case, the positions 
of the particles are more naturally labelled by their angles 8i 's (rather than their positions xi 's) , in 
terms of which the Hamiltonian Hl (|3"9"| associated to this diffusion reads (see |19| ) 

9 2 N pp. 

i=l 1 

which actually comes from the expression of the bi-dimensional Laplacian in terms of polar variables 
(we recall that the radius of the circle is L/2ir and the diffusion coefficient is D = 1/2). It is easy 
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to see that the one particle eigenfunctions of (|53p 4> n {0) which satisfy the boundary conditions 
4> n (0) = 4> n {B + 2tt) are given by 

if^^^eN, (54) 

V Z7T 

such that the N particle anti-symmetric eigenfunctions \Pe{6) and associated eigenvalues E are 
given by [see Eq. (|4*5|) ] 

1 2tt 2 N 

* e{6) = 7m. i<t<N ^ ^ ' e = jtI1 n * ■ ( 55 ) 

* ' i—l 

After some manipulations similar to the one performed before in the case of model A, one arrives 
straightforwardly at the formula for Gn{L) given in Eq. (fl9|) . 



3.3 Model C: reflecting boundary conditions at x = and x = L 

In this case we consider N non-intersecting Brownian motions with diffusion coefficient D = 1/2 
on the segment [0, L] with reflecting boundary conditions both at x = and x = L. The analysis is 
then exactly the same as the one performed for model A in Eqs. (|5§|) - (|4"3"|) except that in that case 
the single particle eigenfunctions 4> n {x) must satisfy: d x (j) n {x) = at x = and x = L. It is then 
easy to see that they are given by 

, . , / 2 /mrx\ , . 

4> n (x) = \l — cos^— — J . (56) 

After some manipulations, exactly similar to the one performed in the case of model A, one obtains 
the formula for En(L) given in Eq. (pH) . 

In the following sections 21 [SJ HO we analyze these formula ([71 [T^l lUTj) in the limit of large N. 



4 Large TV analysis of the maximal height of AT non-intersecting excursions 

In this section, we focus on the distribution of the maximal height of N non-intersecting excursions 
F N (L) given by 

. oo oo N 

ni— — oo tin — — oo i—l 

where Z\jv(2/i, • • • ,Un) is the Vandermonde determinant 

l<j<k<N 

and the amplitude is given by 

n 2N 2 +N 

An = 2- + fn- 1 r ( 2 + j) r(l +J ) ' (59) 
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We introduce the parameter a defined as 



a = 21? ' (60) 



such that F?j{L) can be written as 



F N (L) 



Q(a,N)= £ 



JV 



(61) 
(62) 



ni— — cjo tin — — co i— 1 



The goal of this section is to provide the large N asymptotic analysis of this formula (pTTl 
Ej2J) , In the large N limit, the leading contribution to Fn(L) comes from L ~ C(v^/V), and hence 
a = tt 2 /(2L 2 ) <C 1. It is thus natural to approximate the multiple sum in f2(a,N) (f6"2")l by a 
multiple integral and then evaluate this multiple integral, in the large N lim it, using a saddle point 
method. This saddle point method yields a non trivial result for L < ^/2N, i.e. for the left tail of 
the distribution, and it will be analyzed in subsection 14. II This calculation was done for YM^ with 
the gauge group G = U(7V) by Douglas and Kazakov [55] and we show below that, after simple 
manipulations, this saddle point analysis in the present case where G = Sp(2iV) can be achieved 
without any additional calculation, by borrowing the results of Ref. [29]. On the other hand, for 
L > y/2N, the saddle point result for £2{a, N) in Eq. (|62|) compensates exactly, to leading order, the 
large N behavior of the prefactor / 2+ f / Hfjo r ( 2 + j) r (% + j) in E Q- dSU- Hence for L > V2N 
the net result of this saddle point analysis is simply F^(L) = 1. The analysis beyond this leading 
(and trivial) order requires the calculation of non-perturbative, instanton-like, contributions to the 
partition function. This was first done by Gross and Matytsin [30] in the case of the gauge group 
G = U(iV), and then by Crescimanno, Naculich and Schnitzer [31] for other gauge groups including 
G = Sp(2iV), using the method of discrete orthogonal polynomials. We describe this powerful 
method in subsection 14.21 which then allows us to analyze the right tail of Fn(L) in subsection 
14.31 Finally, this method of orthogonal polynomials also allows us to analyze the central part of 
Fn(L), which describes the typical fluctuations of the maximal height, by studying this system of 
orthogonal polynomials in the double scaling limit, which is done in subsection 14.41 



4.1 Coulomb gas analysis for large N and the left tail 

Physically, it is natural to consider the expression of 0(a,N) in Eq. (|6"2"|) as the partition function 
of a discrete Coulomb gas. The large N analysis can then be performed along the line of the work 
of Douglas and Kazakov in the context of the partition function of the YM2 with the U(A r ) gauge 
group [29], using a constrained saddle point analysis. As we show below, although the discrete 
partition sum f2(a, N) in our model is not exactly the same as that of Douglas and Kazakov, it is 
nevertheless possible to transform directly their result to our case without any further additional 
work. This then allows us to obtain the left tail of the distribution Fj^(L) in the large N limit. 
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For this purpose we first rewrite F^{L) in a slightly different way 

2 N 



F N {L) = 
Q{a,N) 



OO OO iV 



N2 + ^n(a,N), 



(63) 
(64) 



«i=0 



TIM— 2 — 1 



where the integers rij are then all positive ri^ > and recall that a = ir 2 /2L 2 . We now regard the 
summand in Eq. (|64[) as a function of the variables 



2TV 



with i = 1, • 



,TV 



(65) 



In the large TV limit the variables {ni/^TV)},-^... jy approximate the coordinates of a continuous 
TV particles system and associated with the particles is a density p(x). Because of the Vandermonde 
determinant these particles experience hard core repulsion and behave as fermions. The configura- 
tion with highest density corresponds to the case where n\ = 1, n% = 2, • • • , njv = TV for which the 
local density is uniform and given by p max = 2, the lattice spacing being 1/(2TV). Therefore the 
density p(x) must satisfy the constraint 



i N r 
p( x ) = j^z2^ x - x ^ <2, y o dxp( x ) = i 



(66) 



where the average (■ ■ • ) is taken with respect to the discrete weight in Eqs. (jBlJlIM)) and where [0, a] 
is thus the support of the mean density, a remaining to be determined. 

The main idea is that in the large TV limit the discrete sum can be replaced by a multiple integral 
over continuous variables and this multiple integral can be viewed as the partition function of a 
Coulomb gas with Xi = rii/2N denoting the position of the i-th charge. The next step is to replace 
this multiple integral over N variables by a functional integral over the coarse grained density p{x) 
a la Dyson [BT1I621IS5] (for a review see Refs. [541164) ). Skipping details, one finds that in terms of 
the scaling variable 



h = 



L 



'2N 



one can approximate to leading order in the large TV limit, 

Vp(x)e- N ' 2s ^ , 



F N (hV2N) 

_2 



S[f>] 



IT' 
1? 



dxx p(x) 



dx 



dx'p(x)p(x')ln\x 2 ~x r2 \ + C 



p(x) dx — 1 



(67) 



(68) 



(69) 



where C is a Lagrange multiplier that enforces the normalization condition of the charge density. 
If we introduce p(x) = p(x)/2, for x > and p(x) = p(—x)/2 for x < it is easy to see that 



S[p] = 2S DK [p] 



dxx 2 p{x) 



dx 



dx'p(x)p(x')\n\x-x'\+C 



p(x)dx — 1 



(70) 
,(71) 
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where C = C/2 is a constant and Sdk[p) is exactly the action studied by Douglas and Kazakov 
[29] with the substitution h 2 = ir 2 /A. 

In the large N limit, the path integral over p in Eq. (|68|) can be evaluated by the saddle point 
method, giving 



Vp(x)e' N2s ^ ~ cxp {-2N 2 S DK [p*}) , (72) 



where p* is such that 



SS DK [p] 



= , (73) 

p=p" 



5p{x) 

which gives an integral equation for the saddle point density 

—rx 2 ~2 p*(x')]n\x-x'\dx' + C = (74) 

2h J-a 

that holds only over the support x G [—a, a] where p*(x) is nonzero. Taking a further derivative 
with respect to x gets rid of the Lagrange multiplier C" and leads to a singular integral equation 
for p*(x) 

£x-f£Vl& = 0, (75) 
2h J-a x — x' 

where 4" stands for the principal part of the integral, together with the constraints on p(x), inherited 
from the ones on p in Eq. (|66p 

p*(x)dx = 1 , p*(x) < 1 ,\/x G [—a, a] ■ (76) 

If one first discards the constraint p*(x) < 1, the solution of the saddle point equation ([75]) is 
simply given by the Wigner semi-circle law |64| 



7T 



= ^V^"* 9 ' (77) 

and thus a = 2h/n. This density p*(x) has its maximum at x = 0, where p*(0) = therefore the 
Wigner semi-circle la w is the solution of th is co nstrained saddle point equation as long as h > 1. 
Noting that h = L/V2N > 1 means L > \/2N, this solution then corresponds to the right tail of 
the distribution Fn(L) (see Fig. [4]). 

What happens when h < 1, or equivalently L < \f2N7 Clearly, the Wigner semi-circle is no 
longer the right solution since it violates the constraint p*(x) < 1. So, there must be another solution 
to the singular integral equation when h < 1. When h approaches 1 from above, the maximum of 
the semi-circle at x = just touches 1. As one decreases h below 1, the density at x = cannot 
obviously increase beyond 1 since its maximum value is 1. What happens instead is like a 'wetting' 
transition, i.e., the density becomes flat and equal to 1 symmetrically around x — over an interval 
[—b,b] and outside the interval it has a nontrivial shape (see the right panel of Fig. @|. In other 
words, the solution reads 

p*(x) = 1 , x G [-b, b] , p*(x) = p(x) , x G [-a, -b] U [a, b] , (78) 
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-2h/n 2h/ir -a -b b a 

h>l h<l 



Fig. 4 Sketch of the density p*(x) for h < 1 (left tail), and h > 1 (right tail). The Douglas-Kazakov 
transition occurs for h = 1, when L = L C (N) = \^2N in model A, where p*(0) = 1. 



where p(x) is iron zero only for x 6 [—a, —6] U [a, b]. Substituting this ansatz in Eq. (|75p provides 
an integral equation for the nontrivial part p{x) of the density 



7T 



2 



-b\ r p(x>) 



, -x - In - 4 -!-—t-dx' = . (79) 

2/i 2 \x + b) ]^ a x-x' K ' 

Note that both a and b have to be self-consistently determined. This integral equation for p can be 
solved explicitly in terms of elliptic functions [29] (here we give the expression from Ref. [65]): 

p(x) = — V(« 2 ~ x 2 )(x 2 - b 2 )n, (~, - ) , xe [-a,-6]U[o J 6] , (80) 
7rax \ x z a J 

where i7i(y, fc) is the elliptic integral of the third kind, which admits the following integral repre- 
sentation 

1 f 1 dz 1 
n 1 (y,k) = - = — 81 

The constants a and & in Eqs (|78[ 180[) can be expressed in terms of standard elliptic integrals (using 
the notations of Ref. |66| ) 



K( S ) = I ' * - , Efe) = /' *£E& . 

Jo \/l - y z z z vl — z z Jo 

Then a and & are determined by the following relations 



1-z 2 



[2E(fc) - (1 - fc 2 )K(fc)] = 1 (84) 



n 2 



= 4K(fc) . (85) 

Reverting back to our original problem, described by the action S[p] in Eq. (|68p . with /5(x) = 
2p(:r), we get that the saddle point is obtained for p{x) = p*{x) where, for h > 1, 
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p*{x) 



, x < 



Ah 2 



TV 
h 2 

, x > 



2h 



2h 

< x < — 

7T 



(86) 



while for h < 1 one has 



p*{x) 



, x < 
2 , < x < b 
4 



V(a 2 -x 2 )(x 2 -b 2 )IIi 
ixax 

, a; > a . 



b 2 b 



b < x < a 



(87) 



These explicit expressions for p*(x) given in Eq. (|77|) for h > 1 and in Eqs. ([751 180|) for h < 1 
can then be plugged into Eq. (f72|) as done in Ref. [29] to obtain 



where 

3 X 

4 ~ 24 2 

while has a more complicated, though explicit, expression in terms of elliptic integrals: 



F-(X) 



■In A" 



(89) 



12 



(1-fc 2 ) 



24 + 96 



(1 



'.,2\2 



X 



(90) 



where a and fc are defined by the relations in Eqs. (|83[) - (|85| l with the substitution, in the last relation 
([55)1 tt 2 /h 2 — > X. This result in Eq. ([55[) yields the expression of the left tail of Fn(L) announced 
in the section 2 in Eq. ([T2"j). 

The result of this Coulomb gas approach in Eq. ([55j) indicates that this quantity Fn(L) exhibits 
a phase transition at h = 1, i.e. at L = L C (N) = \ 2N. To investigate the order of this transition, 
we need to analyze the behavior of the rate function <p~^(h) in Eq. (|88)1 when h — > 1, which requires 
the behavior of F-(X) — F + (X) for X close to 7r 2 , which was computed in Ref. [55] asQ 



- F + (X) = JL (X - tt 2 ) 3 + 0((X n 2 )*) 



(91) 



This cubic behavior implies a third order phase transition as the third derivative of the rate function 
presents a discontinuity when h crosses the critical value h = 1. Note also that, when h — > 1, h < 1, 
one obtains that F^{hy/2N) behaves like 



F N (hV2N) ~ e" Ar2 ^ (1 -' l)3 , 
which yields the asymptotic behavior of the rate function 4>~^(h) given in Eq. |T 



(92) 



Note that we have corrected a sign error from Ref. |29] and reproduced also in Ref. |19] 
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We see in Eq. (|55|) that this saddle point approach does not say anything relevant about the 
right tail of the distribution. It just gives the trivial leading order result that Fn(L) ~ 1 for 
L > s/2N. To investigate this regime beyond this trivial leading order, we have to find some other 
method. Fortunately such a method exists and it involves discrete orthogonal polynomials which 
were introduced and studied by Gross and Matytsin [30] in the context of YM2 with the U(-ZV) 
gauge group. We can adapt their method to our problem as shown in the next subsection. We will 
see that this method also allows us to study not just the extreme right tail but also the typical 
fluctuations, for L ~ \/2N, which is described by a double scaling limit which we analyze below. 



4.2 General expression in terms of discrete orthogonal polynomials 

To analyze (|6T1 162|) in the large N limit, we follow Ref. [30] and introduce the following discrete 
orthogonal polynomials 

00 

Pk(n\a) Pk ,(n\a)e- an2 = 5 k , k .h k {a) , (93) 

TL— — 00 

where we recall that a = ir 2 /2L 2 . Here p k (n\a)'s are monic polynomials of degree k: 

p k {n\a) = n k + ■ ■ ■ (94) 

and hk((x)'s are to be determined. 

These orthogonal polynomials turn out to be very useful as they allow us to express F^(L) in 
a rather compact form, as we show now. First notice that 

n 

ln i A(n 2 ,n 2 ,,---,n 2 N )= det Jpa-i^j)] • ( 95 ) 

Ki,i<N 



Using this in Eq. ([p^|) we get 



00 

e" 

n 1 — — OO 7lN = — QO 



n(a,N)= £ ... £ 



, det r [P2t-i(%)] 

l<t,3<N 



(96) 



Next, using the definition in Eq. (|93|) and the discrete version of the Cauchy-Binet formula one gets 

N 

n(a,N) = r(N + l)Y[h2j-i(a). (97) 
3=1 

Thus, if we can determine the /ifc(a)'s, in principle we can determine f2(a,N) exactly. 

So, our next task is to determine the /ife(a)'s. As we show below, one can actually set up a 
very nice recursive procedure to determine the /ife(a)'s. To proceed, we notice that the two first 
polynomials are obtained straightforwardly as 

p (n\a) = 1 , p\{n\a) = n , (98) 

from which one gets the first amplitudes 

00 00 
h (a) = J2 e ~ an2 ' h ^ a ) = n 2 e~ an2 . (99) 
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Note that given that the weight e an is symmetric under the transformation n — > —n and given 
that the sum in Eq. (|93|) runs over neZ, one can show (for instance by induction over k) that 

P2k (n) = n 2k + a 2k -2n 2k - 2 + ■ ■ ■ + a a (100) 
P2k+i(n) = n 2k+1 + b 2k -in 2k - 1 + ■ ■ ■ + bin . 

Because p k 's are orthogonal polynomials (j9"3")l they satisfy a three term recursion relation [6TJ : 

np k (n\a) = p k+ i(n\a) + S k {a)p k (n\a) + R k (a)p k -i{n\a) . (101) 

Multiplying the above recursion relation (jlOip by p k (n\a)e~ an and summing over n, one obtains 
from ([TOO"]) that 

S k (a)=Q. (102) 

On the other hand, multiplying both sides of Eq. (|101[) by p k ^i(n\a)e~ an and summing over n 
one obtains 

R k{a ) = . (103) 



From Eq. (|99|) . one has 

fil(a) = Ma) = "e 55 — = "^ ln/lo(a) ■ (104) 

Inspired by the relation found above for i?i(a) (|104p . we differentiate the orthogonality condi- 
tion (|93|) with respect to a 

-^-h k (a) = ~(np k \np k ) + 2{p k \d a p k ) , (105) 
aa 

where we have introduced the notation (f\g), for two functions /, g: 

oo 

(/Iff) = E /(")5(^)^ Q " 2 ■ (106) 

71= — oo 

Using the fact that p k is a monic polynomial (fMf . <9 Q pfc is thus a polynomial of degree at most 
fc — 2 (jlOOp . And therefore, from the orthogonality condition ([93| . the second term in (| 105[) vanishes. 
Finally, using the recursion relation (|101j) one obtains 

^h k {a) = - (R k+l (a)h k {a) + R k {a)h k {a)) . (107) 
aa 

From the relation found above for Ri(a) (|104p we see, from Eq. (|107p that one has, for consistency: 

R (a) = . (108) 
By subtracting the relations f| 107[) for -j^h k (a) and -^h k ^x(a) one obtains 

- -^-lni? fe (a) = R k+1 {a) - i? fc _i(a) . (109) 
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Let us then summarize the main program to be carried out. We have a set of variables i?fc(a)'s 
that satisfy a Volterra type recursion relation in Eq. (|109j) . starting from the initial values: Ro(a) = 
and Ri(a) = — -4-ln/io(a0 where ho(a) = X^^L-oo e ~ Q ™ • This is a closed set of deterministic 
recursive relations, albeit nonlinear. In principle, if we have the solution Rk(a) for arbitrary fc, then 
we can determine the ft,fc(a)'s by iterating the relation () 1 03[) : hk(a) = Rk{ot)hk-i{a), starting from 
ho(a) = 5^^L-oo e_Q ™ ■ Once we have the /ife(a)'s for all k, by taking their product in Eq. (|9"T)) . we 
can determine, at least in principle, f2(a,N) and hence subsequently Fn(L) using Eq. (|6ip. 

Actually, to get rid of the complicated prefactor in Eq. (pTj) . it turns out to be useful to consider 
the following ratio 

F N+1 {L)F N -i(L) a 2 h 2N+1 {a) a 2 

fJiL) = N(N+ 1/2) h 2N ^(a) = NWTW) R2N+lia)RMa) ' (110) 

where we have used Eqs. (|61"]) . (|97|) and t]103[) to simplify. This will turn out to be the most useful 
form on which one can perform the asymptotic analysis for large N and large L. 

So, essentially, our aim and the main challenge are to solve the nonlinear recursion relation for 
the i?fc(a)'s in Eq. (|109|) . Evidently, this can not be solved for all k exactly. However, for large 
k and large L (i.e., small a — ir 2 /2L 2 ), it turns out that one can make some scaling ansatz and 
solve the resulting differential equation for the scaling functions. This provides us a handle on the 
asymptotic analysis of Fn{L), as we show below in the next two subsections. 



4.3 Large deviation regime: right tail 

2 

To begin with, we first consider the limit when L — > oo, i.e. a = yp- — > 0, with N fixed to a large 
value. To get insight into this limit, it is useful to study the amplitudes ho(a), hi(a) in Eq. (|99|) 
when a — > 0. For this purpose, we use the Poisson summation formula: if for a given function fix) 
we define its Fourier transform 

f(q) = ^J° f(x)e^dx, (111) 

then one has 

oo oo 

E /(»)= E /( 2 ™)- ( 112 ) 

7i— — oo n— — oo 

When applied to h (a) and hi (a) this Poisson formula yields 

W = vf E e ^ = \/I (l + 2e-£+0(e-^)) , 

71— — OO 

hl (a) = ^ f e -^(-2n 2 n 2 +a)=A(l-(4^-2)e-4 + o( e -^ 



(113) 



n— — oo 

from which one gets 



Ri(a) = - ^e"^ + 0{e~^/a) . (114) 



2o 



More generally, one can show that, to leading order, Rk(a) is given by 

fl*(«)-^> ( 115 ) 

which indeed satisfies the above Eq. (|109[> . Motivated by the result for R\{a) (|114j) we assume that, 
when a — > 0, the leading correction to (|115p is of the form 



^(«) = T^+c fc (a)e-^+0(e- 1 S i /a), (116) 
la 

such that Ck(a)e~~ <C 1. Note that from (|108p and Q114p one has 

co(a =0, ci (a) = -2-= . (117) 

ar 

By inserting this ansatz (|116p into the above recursion relation satisfied by Rk(a) (|109[) . assuming 
that Cfc(a)e _ ~ <C 1, one finds the following recursion relation for Ck(a): 



2 2a 2a n 2 

c fc+ i(a) - Cfc_i(a) = -r c fe( a ) ~ "i" c fc( a ) ~ y^ c *( a ) • ( 118 ) 



Following Ref. [3D] we introduce the new variable 



£ = - (H9) 



and rewrite 

Cfc(a) = - ^ 2 G fe (0 , (120) 
such that Gfc(£) satisfies the recurrence relation, obtained from (|118[) . 

G k+1 (0 - G fc _x(0 = H [£G' fc (0 + G fc (0 - • (121) 

From Eqs. (|117p and (|12ip one obtains explicit expressions for the first Gfc's as 

G (0 = , Gi(0 = 1 , G 2 (0 = -2^ + 2 , G 3 (0 = 2£ 2 - 6£ + 3 . (122) 
This recursion relation (|121|) can be solved by introducing the generating function 

oo 

G&z) = J2Gk(0z k , (123) 
fe=i 

which satisfies the partial differential equation (pde) 

<• -«*>!; -«f-H-'Cf-u)e. (124) 
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where we have used the short notation G = G(£,z). This pde (jl24[) admits the following exact 
solution [30] 



z ) = 77— TI CX P ( "A ) , (125) 



2£z 

(1 - z) 2 V ~ z 
from which Gfc(£) is readily obtained as 

dz 1 z.,^ ^ f dt __ 2 £t i - 1 



where we have performed the change of variable £ = z/(l — z). In the above expression the contour 
C z encircles the origin z = in such a way that z = 1 is outside of the contour while the contour 
C* encircles the origin t = and passes to the right of t = — 1. One recognizes from Eq. (|126|) . for 
instance by expanding (1 + l/t) k in powers of 1/t and integrating term by term that Gfc(£) can be 
written as 



fc-i 



= E 7 mw/ 7TT(20 fc = 4-i( 2 , (127) 

z — ' m [m + 1 ! /e — m — 1 ! 

m=0 v ' v ; 

where L^^x) are generalized Laguerre polynomials |67| . Therefore one gets, coming back to the 
variable a = (|119|) . when a — s- 



i? fc (a) = — - 2^LW — + G(e- 4 ^/«) . (128) 



2a "V V " 

We now want to analyze this formula for i?fc(a) (|128p when k ~ 0{N) is large, see Eq. (jllOp . 
For this purpose, we analyze the integral formula for Gfc(£) in Eq. (|126p using the saddle point 
method and write 

G fc (0= / A e -«H-kln(l+l/t) (129) 
/c t 2tti 

which we analyze in the limit where k, £ 3> 1 and iV ^> 1 such that 

k £ 2L 2 . . 

with y and T kept fixed. Therefore from Eq. (|129[) one has 

<w=£ ^ e_Ars(t) ' 5 w = 2T< -^ ln ( 1+ i) ■ ( 131 ) 

which can then be evaluated by a saddle point method when N — > oo. This saddle point is such 
that dS/dt\t=t* = 0, where S'(i) reaches its maximum on the real axis, yielding 



e = + . (132) 



It corresponds indeed to a maximum of S(t) on the real line 



If* 4- 1 

S"(t*) = -y . . < . (133) 
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Note that this saddle point solution (|132[) makes sense if and only if t* is real, which implies 



2 |<l^fc<| 



(134) 



Now given that the largest value of k involved in the expression of F^ ( N) in Eq. piOp is k = 2 TV + 1 , 
this means that this saddle point analysis holds for 



L > V2N , 

which thus allows to study the right tail of the distribution. This yields finally 



G k (0 



1 



-NS(f 



2tt y/N\S"(r)\ 



(135) 



(136) 



where the minus sign in Eq. (| 136[) comes from the fact that the integral over the contour Ct in Eq. 
(11291) is counterclockwise oriented, which thus runs, close to t* from 3(t) > to 3(t) < and hence 
the minus sign. One finally obtains, plugging the explicit expression of t* (|132p into Eq. (I13(ip 



G k (Oe 



\k-l 



cxp 



(2k 



y(x) = \/l-x - - In 



1 + s/T^X 



1 - jT^x 



(137) 
(138) 



where the prefactor (— l) k in Eq. (|137[) comes the fact that 1 + 1/t* < such that when one evaluates 
e ~NS(t ) jjjgj.g j s a (— l) fc prefactor coming from exp[fcln(l + 1/t*)] (|13ip . Note in particular the 
asymptotic behavior of j(x) 



7(«) ~ oC 1 



2; — >• 1 



(139) 



which will be useful in the following to understand the matching between large and typical fluctu- 
ations. Finally one has in this regime (|130[) : 



R k (a) = ^ + {-l) k 



2 1 /2 7r 5/2 



2k 

T 



-1/4 



exp 



u 



(140) 



where we recall that £ = it 2 /a = 21? . We now want to use this asymptotic behavior (|140[) to study 
the right tail of the cumulative distribution Fn(L). For this purpose, we start from the expression 
given in Eq. (|?7|) to write (see Appendix Al) 



F N (L) 



r(N + i) 



N-l 



+ 3 



■a NL+ f [hoWR^a)]" ] [ [R 2k {a)R 2k+1 (a)[ 



N-k 



(141) 



fe=i 



with R k (a) given in Eq. (TTT6| . which for k ~ O(iV) reduces to Eq. ([T3D]) . To study this ex- 
pression ()141[) in the large deviation regime, we write R k (a) as in Eq. f| 1 and treat the term 

2 

c k (a)e~~ in perturbation theory. After cumbersome manipulations left in Appendix |A] we arrive 
at the following formula 



In F N (L) 



,-Z 



[G 2 n(0+I2n(0] + O exp 



-2^7 



2N 

T 



(142) 
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in terms of the variable £ given in Eq. (|119p and the functions G2at(£) given in Eq. (|127p and I2n(Q 
given by 



N-l 



G 



2k+l 



(0 



2/c + 1 



(143) 



This formula in Eq. (|142j) yields back the formula obtained in Ref. [31] [see their Eq. (28)]. Us- 
ing the integral representation for Gk{£) obtained before in Eq. (|126[) one can obtain an integral 
representation for l2iv(£)> see Appendix IX| under the form 



2H 



C t 



dt e-« 
2ttI 1 + It 



1 

1 + 7 



2A ! 



(144) 



The analysis of I2n(Q in the large deviation regime (|130|1 can then be carried out as before for 
G2n{0- By comparing the two formulas (|126p and (| 144|) it is easy to see that in the large L and 
large N limit one has 



1 



l + 2f 



-G 2JV (£) = 



1 



V 1 - 4 / T 



G2n(C) 7 



(145) 



where we have used the expression for the saddle point in Eq. (|132jl . Finally one obtains that in 
the large deviation regime 



1 



1 



v/1 - 4iV/C 



1 



2^ e 



i - 



4iV 



-1/4 



exp 



(146) 



which is negative, as it should, since Fn(L) < 1. Recalling that £ = 2L 2 one obtains from Eq. (|14G|) 
the right tail behavior announced above in Eq. (fl~2|). 



1 - F N (L) ~ cxp -N(j)+(L/V2N) 
(f)+(x) = 4x 2 7(1/.t 2 ) = 4.ta/x 2 - 1 - 2 In 2x (V2; 2 - 1 + ^ 



- 1 



(147) 



where we have used the explicit expression of j(x) in Eq. (|137p . From the asymptotic behavior (|139p 
one gets immediately 



+ (z) = ^(z-l) 3 / 2 + o((x-l) 5 / 2 ) ,x^V 



(148) 



which yields the asymptotic behavior announced above in Eq. (|13|) . Note in particular that if we 
set L close to \/2N such that 



L = V2N + cN-*s , c= 2~ 11/6 
one has, using this asymptotic behavior in Eq. (|139[) 



lnFjvOL) 



1 



4v^s 3 / 4 



_2,3/ 2 

e 3 



iV-i/3 



1 



e -§« 3/2 +0 (7V- 2 / 3 ) 



(149) 



(150) 



which is useful to study the matching, discussed in section 2, between the right tail and the central 
part of the distribution, which describes the typical fluctuations around L = \J2N. 
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4.4 Double scaling regime 

We now focus on the typical fluctuations of the top path and study Fn(L) for L close to y/2N. 
This corresponds to a double scaling limit where both L and N are large but such that L — \f2N = 
0{N~ 1 ^). Hence we set 

L = V2N + cN- 1/6 s , c = 2- 11 ? 6 . (151) 
We define a "running variable" x k as 

x k = n% 3 ( 1 - — ) , n cr = £ = L 2 . (152) 
\ n cr J 2 

The goal is to analyze the recursion relations for R k {a) in Eq. (|109[) when fc is also of order O(N). 
This can be analyzed in this double scaling limit by assuming the following ansatz [SU] (which was 
later on proved rigorously in Ref. |46| ) 

2 

R k (a) = % + n% 3 f+{x k ) + n% 3 f+{x k ) + n cr f+(x k ) + 0{n% 3 ) , k even , (153) 

7T 

2 

R k (a) = + r# 3 /r(a*) + ni /3 / 2 ~(^) + n cr f^(x k ) + 0{n 2 J 3 ) , fc odd . 

The fact that we assume different functions , , depending on the parity of k is guided by 
the previous result ()140[) where we showed that the leading correction comes with an amplitude 
proportional to (— l) fc , hence depending on the parity of fc. By substituting this ansatz (|153|) into 
the recursion relation in Eq. (| 109)) one obtains that 

= -/f(z) = ~h(x) , (154) 
where the function f% satisfies a Painleve II equation (PII) 



o5/3 

f['(x) = Axhix) + -ff(x) , h{x) r Ai(2 2 / 3 x) , (155) 



where A\(x) is the standard Airy function. Note the minus sign in the asymptotic behavior which 
is missing in Ref. [3U]. It can be expressed in terms of the Hastings-McLeod solution of PII q{s) in 
Eq. © as 

h{x) = -?-l q (2 2 / 3 x) , (156) 
n z 

where 

q "( S )=2q 3 ( S ) + sq (s) 1 q(s) ~ Ai( S ) ~ _L^-fs 3/ \ (157) 
One can further show [3D] that the functions , / 2 ~ i n Eq. (|153|) satisfy 

+ / 2 -(x) = + ~yfi( x ) . (158) 

7T Z 2 
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and that the functions /g , / 3 in Eq. (|153[) satisfy 

f+(x) + f 3 -(x) = ixA(x) + i 2 /iW/ 2 1x) - y A 3 (x) + . (159) 

Using this ansatz (|153p . one can then analyze Fn(L) in the double scaling limit. For this purpose, 
it is useful to start from Eq. (jllOp : taking the logarithm on both sides and performing the large n cr 
expansion of the right hand side, one obtains after a tedious though straightforward calculation, 
using (fT54 l fT58 | : 

hxF{L,N+ 1) +lnF N _ 1 (L) - 2lnF N (L) = -n^' 3 ^- (f?(x 2N ) + ^fi(x 2N )) (160) 



On the other hand, if we assume the following ansatz for hi F^(L) 

]nF N (L) = Y(x 2N ) + n- r 1 / 3 H(x 2N ) + 0(n~ 2/3 ) , (161) 
the expansion of the left hand side of Eq. (|160|) yields 

\nF N+1 {L) + lnFjv-i(L) -2lnF N (L) = 4, n - 2 ' 3 Y" {x 2N ) + 4n^H"{x 2N ) + Ofc 4/3 ) ■ (162) 
Therefore identifying the terms with the same power of n cr in Eqs. (| 160[) and (|162p one obtains 

4Y"(x) = ~ (h(x) 2 + J>/i(z)) (163) 

m"{x) = ^-( K f 1 {x)f[{x) + ^fl{x)\ . 

In terms of the variable s defined in Eq. (|151|) one has from (|152[) 

x 2N = 2- 2/3 s , (164) 

and if we define 

Y{x) = Y(2 2/3 x) , H(x) = H(2 2/3 x) , (165) 

the above equations (|163[) can then be written in terms of q(s), the Hastings-McLeod solution of 
PII © as 

Y"(s) = -±( q 2 (s)- q i( S )) (166) 
H"(s) = 2-V3 Us) q '(s) - = -2- 1 /3 y - (s) . 

Therefore, from Eqs (|161[) . (|165p together with Eq. (|166p one obtains, using that Fn(L oo) = 1 



\nF N {V2N + 2- 11/6 N- 1/6 s) = Y(s) - N~ l/3 2- 2/3 Y' (s) + (D(N~ 2 / 3 ) (167) 

2 



1 f°° 

Y(s) = -- (x-s)q 2 (x)dx-- I q(x)dx. (168) 
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So that finally one has 

F N (V2N + 2~ 11 / 6 N~ 1 / (i s) = J"i(s) - N-V 3 2- 2 / 3 T[(s) +0{N~ 2 ^) . (169) 

Note that using the large s behavior of q(s) in Eq. (|157j) one obtains the large s behavior of .Fi(s) 
as: 

lnJ-x( S )~--y 9 ( s )d a: ~- 5 -^ 57 j e -t' 3/a (l + 0( s - 3 / a )) J (170) 
and therefore, expanding each term of the large N expansion in Eq. f|169[) for large s, one gets 

\nF N {V2N + 2- u ^N~ 1 ^s) = - - * — e^'fl + 0( S ~ 3 / 2 )) (171) 

- ^" 1/3 2873^174 53/2 (! + °(^ 3/2 )) + (^" /3 ) ' 

This expansion matches perfectly with the large deviation behavior obtained above in Eq. (|150p . 
We notice that this is not only true for the leading term of order O(N ) but also for the subleading 
one, of order 0(N^ 1 ^ 3 ). 

Although the leading term in Eq. (|169p was already obtained in Ref. [T!5] , and subsequently shown 
rigorously in |3B], we also obtain here the first correction, of order 0{N~ X / 3 ) to T\. Note that this 
correction, being proportional to J\{s), can also be written as a simple shift of the argument of T\ 
as 

F N (V2N + 2- 11/6 N- 1/e s) = ^{s ~ N- 1/3 2~ 2/3 ) + (D{N- 2/3 ) , (172) 

such that this leading correction only affects the first moment, the finite N corrections to the higher 
cumulants being at least of order 0(N~ 2 / 3 ) (see also Ref. [68] for a related computation in a sightly 
different context). 



5 Large AT analysis of the reunion probability of TV non-intersecting Brownian 
motions on a circle 

In this section we focus on the ratio between reunion probabilities Gn(L) for non- intersecting 
Brownian motions on a circle, i.e. with periodic boundary conditions. We start with the expression 
given in section 2 in Eq. (|19jl 

D oo oo 

Gn(L) = ^ £ ... £ A%(n u ...,n N )e- 2 -£^n> (m) 

ni- — oo n^ — — oo 

where An{jji, • • • , Un) is the Vandermonde determinant (|58[) and the prefactor Bn is given by 

1 

(27r)^-^n^r(i + 2) 
Here we introduce the parameter a defined as 

2?r 2 

a = — , (175) 



B * = 7Z — n ^W^T. . ^ ■ (174) 
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and rewrite 

Gn(L) = — ^ — = V ••• V 2\jv(ni,--. ,njv) 2 e- a ^f=i"' . (176) 

Note that, for simplicity, we use the same notation a as in model A, and in model C below, although 
its relation to L differs from one model (|f)0"|) . (|199|) to another (|175p . 



5.1 Coulomb gas analysis for large N and the left tail of Gn{L) 

This case corresponds exactly, up to a multiplicative prefactor, to the partition function of Yang- 
Mills theory on the sphere with the gauge group U(iV), which was analyzed by Douglas and Kazakov 
in Ref. [29]. The transposition of their results to the present study was already presented in Ref. 
|19| . There it was shown that the transition for Gn(L) happens for L = 2yN such that, if one 
introduces the parameter r = Lj2\[~N one has, similarly to what we have obtained before for model 
A " 



1 , r > 1 

lim -— lnG JV (v / 4iVr) = <^ ' " , , . , (177) 



where the function <f> A is defined in Eq. (|88| . Here, the asymptotic behavior of Gn(^VN t) for r 
close to 1 is thus 

G N (VlNr) ~ exp f-^N 2 (l - r)») . (178) 



These results in Eq. (|177|) and Eq. (|178p yield the behaviors announced in section 2 in Eqs 
(|2"3")l . for the left tail of Gn{L). However, as before for model A, this saddle point method does not 
give anything meaningful for the right tail, which can be analyzed using the method of discrete 
orthogonal polynomials described before. 



5.2 Analysis of the right tail of Gn{L) 



The right hand side of Eq. (|176|) can be conveniently expressed in terms of the orthogonal poly- 
nomials pk(n) introduced and studied above (|93|) where a is now defined by (| 1 75[) . Using standard 
manipulations (in particular the Cauchy-Binet formula), one arrives at 



JV-l 



G N (L) = Ni aN /2/ ^ + 1 > n Ha) 



(179) 



As we have done before (|141j) . this product of the amplitudes hj in Eq. (|179p can be re- written in 
terms of the Rk's. Substituting then the expression for R^s in Eq. (|103j) and treating the terms 

Ck{a)e~~ in perturbation theory, one obtains (see also Ref. |30| ) 



\nG N {L) = 2e-tG N {0+O ( cxp 



-2^7 



2N 

T 



(180) 
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in terms of the variable £ = Tr 2 /a = L 2 /2, reminding that a = 2ir 2 /L 2 in this case (j!75[) . This 
yields, using the expression for Gjv(£) given in Eq. (|13T() : 



lnGjv(L) 



(-1) 



N-l 



1 



2N 




2 v<V £ 
(_i)*-i 



2iV 



-1/4 



2ttL 




1 



AN 

T 2 



-1/4 



exp 



exp 
■ L 2 



2 7 V L 2 



(2N 

AN 



(181) 
(182) 



which shows an interesting oscillatory behavior with N: this fact is not problematic as Gn(L) does 
not have the meaning of a cumulative distribution. One thus obtains from Eq. (|181|) : 



G N (L) = (-If exp -N0+(L/VAN) , <j>+{x) = 2x'j 



(183) 



as announced earlier in section 2 (|22l) , where we have used the expression of 4>\ (x) given in Eq. (|147l) . 
The asymptotic behavior of cj)^ (x) when x — > 1 follows immediately from Eq. (| 148[) , as announced 
in Eq. (f2"3")l . Therefore, if we study the regime of L close to 2\/]V and set 

(184) 



L = 2^N + 2- 2/3 N- 1/6 t , 
one obtains from Eq. (j!81|) 



1 



2 4 /3 0rfl/4 



_2.3/2 

e a 1 



(185) 



where we notice that this term is of order 0(N^ 1 ^ 3 ): this is responsible for the (mis-) matching with 
the central part of the distribution, which is instead of order O(N ), and described by a double 
scaling limit which we now focus on. 



5.3 Double scaling regime: analysis of the central part of the distribution 

To analyze the central part of Gn(L) we start with the following identity, analogous to the one in 
Eq. (jllOp for model A, which reads here 



G N+1 (L)G N ^(L) 2a h N 2a _ . , 

= = = — Hn(oi) 

[G N {L)] 2 N Iin-! N Ny ' 



(186) 



which is very useful for asymptotic analysis. Using the asymptotic expansion of Rpf(a) in Eq. (|153p . 
one can expand the right hand side of Eq. t]186[) as 



where n cr is given here by 



I? 
4 



(187) 



(188) 
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On the other hand, guided by the result which we have obtained for the right tail of Gn(L) and 
also by the expansion above (|187p . one assumes the following ansatz for In Gn (L) in the double 
scaling regime 

In Gn(L) = Y(x N ) + {-l) N - x n-V*H(x N ) , (189) 
such that the left hand side of Eq. (|186[) admits the following expansion 

In Gjv+i (L) + In Gn-i(L) -2\nG N (L) = (-l) N n- 1/3 4H(x N ) + n - 2/3 Y"(x N ) + ©(n" 1 ) . (190) 
Therefore, identifying the different powers of n cr in Eqs. f| 187[) and (| 190)) one finds 

Y"(x) = -^h(x) 2 , H(x) = ~h{x) . (191) 

If we set 

Y(x) = Y(2 2 / 3 x) , H(x) = H{2 2 / 3 x) , (192) 

F"( S ) = -q 2 (s) , H(s)=2-V 3 q(s). (193) 

One thus obtains 

\nG N (2VN + 2- 2 ' 3 sN- 1 / 6 ) = Y(s) + (-^^"^-^^"^^(s) + 0{N- 2 ' 3 ) (194) 

Y(s) = - {x- s)q 2 {x)dx , (195) 



one has from Eq. (|191[) 



which yields 



lim G N (2VN + 2- 2/3 sN- 1/6 ) = F 2 {s) . (196) 



The result in Eq. (|190[) yields the formula given in Eq. (|27|) in section 2. At variance with model A 
(and model C discussed below) the competition between the two terms ^(s) ~ exp(— |s 3 / 2 ) and 
q(s) ~ cxp (— |s 3 / 2 ) yields an unconventional crossover regime discussed in section 2 and illustrated 
in Fig. [3l which is peculiar to this vicious walker problem. 



6 Large N analysis of the reunion probability of N non-intersecting Brownian 
motions with reflecting boundary conditions 

In this section we focus on the ratio of reunion probabilities En(L) for non- intersecting Brownian 
motions on the segment [0,L] with reflecting boundary conditions at x — and x = L. We start 
with the expression given in section 2 in Eq. (|31l) 



„ oo oo 

E N (L) = I ^ ZW E ■•• E ^{n\,...,n 2 N )e-^U^ , (197) 
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where An(vi, • • • , Un) is the Vandermonde determinant (|58[) and the prefactor Cn given by 

2N 2 -N nN/2-N 2 

Cn = — m — ? , (198) 

nf=" 1 r{2 + j)r(i/2+j) 

ensures that En(L — > oo) = 1. Here we introduce the parameter a given, as in model A (|60p . by 



a = ^ > ( 199 ) 



and write 



E N (L) = N . ° V ... V ^ 2 (n;,-»n^)e- aE "= in ". (200) 



-oo njv — — oo 



6.1 Coulomb gas analysis for large TV and the left tail of En(L) 

The Coulomb gas analysis of En{L), which amounts to the study of the path- integral (over the 
density p) entering the expression of E^{L) as in Eq. (|68p . is exactly the same as the one done 
for Fisf(L). Indeed the difference between Fj^(L) and En(L) is the presence, in the expression for 
Fjy(L), of the product YiiLi n 1- But this term does not contribute to S[p] as hi nf ~ O(N), 
and are thus subdominant compared to 0(N 2 ) terms which contribute to S[p\. Therefore one gets 
immediately the result for the left tail of E^(L) 

lim --L\nE N {hV2N) = i ° ' /? N ~ 1 , „n , ( 201 ) 

where the rate function 4>^(h) is given in Eq. (|88|) . This yields the result for the left tail of En(L) 
announced in Ref. (l34l). 



6.2 Analysis of the right tail of En(L) 

The expression above ()200[) can be conveniently expressed in terms of the orthogonal polynomials 
Pk(n) introduced and studied above (|93|) where a = ir 2 /2L 2 (|199p . Using the identity 

A N (nl ■■■ ,n%) = det [pw-afo-)] , (202) 

l<2J<iV 

and the Cauchy-Binet formula, one arrives at 

E N (L) = TT ky- 2 (a) ■ (203) 

nf=i n2+i)nv2+i) J i i 

As we have done before (|14ip . this product of the amplitudes hj in Eq. (|203p can be re- written in 
terms of the i?fc's. Substituting then the expression for Rk's in Eq. (|103|) and treating the terms 

2 

Cfe(a)e - ^ in perturbation theory, one obtains [see also Ref. |31] ] 



lnE N (L) = e-« [Gjv(0 ~ h N (0) + O exp 



-2£ 7 



27V 

T 



(204) 
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in terms of the variable £ = tt 2 /a = 2L 2 , where we have used a = tt 2 /2L 2 (|199p . Note that this 
expression is very similar to, albeit different from, the expression obtained for lnFjy(L) in Eq. (11421) 
(the term I 2 n{£,) comes here with a minus sign). One thus gets, from the previous analysis (|14GD , 



lnE N (L) 

Therefore, one obtains from Eq. (|205[) 



1 VN / 4AT 



-1/4 



cxp 



<1 



AN 



1 - E N {L) = - cxp [-N</>+(L/V2N)] , 4,+ (x) = 4x 2 7 -j ) = 4>\(x) 



(205) 



(206) 



as announced in Eq. (|34|) . The asymptotic behavior of ^ J (x) when x — > 1 follows straightforwardly 
from Eq. (|148p . Note in particular that if we set L close to y/2N such that 



L = VM + cN'h , c = 2~ 11/6 , 
one has, from Eq. (|205[) . using the asymptotic behavior for 7(x) in Eq. (I139[) 



In E N {L) 



1 



2 3/2 -i /q 



1 



e -i« 3/2 + 0(Ar-2/3) 



4^Fs 3 / 4 2 8 /30Fs 1 /4 
which is useful to study the matching with the typical fluctuations around L = y/2N. 



(207) 



(208) 



6.3 Double scaling regime: analysis of the central part of the distribution 

To analyze the central part of En(L) we start with the following identity, analogous to the one in 
Eq. (|110p for model A, which reads here 



En+i(L)En-i{L) 
[En(L)} 2 



a 



N(N - 1/2) h 2N _ 2 N(N- 1/2) 



R 2N (a)R 2N -i(a) 



(209) 



which is very useful for asymptotic analysis. Using the asymptotic expansion of Rn(ck) in Eq. (|153p . 
one can expand the right hand side of Eq. (|209[) as 



In 



a 



N(N- 1/2) 



R 2N (a)R 2N - 1 (a)j = -n cr 2/3 y (^f 2 (x 2N ) - -^f{(x2N)j 



(210) 



+ n crY (-h(x2N)f[(x 2N ) + ^/l(S2AT)) + 0{n~J' S ) 



where, here 



L 



On the other hand, if we assume the following ansatz for \tiEn(L) 

]nE N (L) = Y(x 2N ) + n- r ^ 3 H(x 2N ) + 0(n; 2 ' 3 ) , 



(211) 



(212) 
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the expansion of the left hand side of Eq. (|210[) reads 

hi E N+1 (L)+\nE N ^{L) -2\nE N {L) = 4n~ 2 / 3 Y"(x 2N ) + 4n~ 1 H"(x 2N ) + 0(n~ r ^ 3 ) . (213) 
Therefore identifying the terms in Eqs. (|210|) and Eq. (|213p with the same power of n cr one obtains 



-y (AO*) 2 - J?/I(* 

iH"(x) = ~(f 1 (x)f{(x)-^f{'(x) 



(214) 



In terms of the variable s defined in Eq. ()207[) one has from (|152|) 
and if we define 

Y(x) = Y{2 2 ' 3 x) ,H(x) = H{2 2 ' 3 x) , 



(215) 



(216) 



the above equations f)2 14[) can then be written in terms of q(s), the Hastings-McLeod solution of 
PII © as 



Y"(s) = --( q 2 (s) + q'(s)) 

H"{s) = 2- 1 / 3 (q(s)q'(s) + lg») = 2-^ 3 Y'"(s) 



(217) 



Therefore, from Eqs. (|212p and (|216|) together with Eq. (|217j) one obtains, using that En(L 
00) = 1 



ln^jv(V2iV + 2- 11/6 7V~ 1/6 s) = F(s) + N~ 1/3 2- 2/3 Y' {s) + 0(N~ 2/3 ) 

lf°° 1 f°° 

Y{s) = --j (x - s)q 2 (x)dx + - q(x)dx . 



(218) 
(219) 



So that finally one has 



E N (V2N + 2- 11 / 6 N- 1 / 6 s) = ^44 + N- 1 '^- 2 / 3 ^- 

Fi{s) as 



T 2 (s) 



Ms) 



+ 0(N" 2/3 ) 



(220) 



where we have used the expression of T 2 and T\ in terms of q(s) given respectively in Eq. ([I]) 
and f4]) . As noticed previously for model A (|172[) , the first correction in Eq. (|220[) , proportional to 
A -1 / 3 , can be absorbed by a shift of the argument s — > s + 2~ 2 / 3 N~ x / 3 . An analogous structure 
is known for the large N expansion at the soft edge of the Laguerre ensemble in random matrix 
theory (see [64], eq. (7.162)). 
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7 Conclusion 

To conclude, we have performed in this paper a systematic study of three different models of N 
non-intersecting Brownian motions on a line segment [0,L] with three different types of boundary 
conditions at x = and x = L: absorbing (model A), periodic (model B) and reflecting (model C) 
boundary conditions. In each of these models we have focused on a normalized reunion probability 
which, in model A, can also be interpreted as the maximal height of N non-intersecting Brownian 
excursions on the unit time interval. We have presented a self-contained derivation of the formulas 
for finite N for these reunion probabilities, whose expressions had been given without the details in 
previous publications in Ref. [T2], for model A, and for model B and C in Ref. [T5]. An interesting 
property of these reunion probabilities is that they are, up to a multiplicative pre-factor, identical 
to the partition function of Yang-Mills theory on the sphere with a gauge group G which is selected 
by the choice of boundary conditions: Sp(2A) for model A, U(JV) for model B and SO (2 A) for 
model C. As a consequence of this correspondence, these reunion probabilities exhibit a third-order 
phase transition, akin to the Douglas-Kazakov transition in YM2, as the size of the system L crosses 
a critical value L C (N) oc y/~N between the left tail L < L C (N) and the right tail L > L C (N). In 
the central part, L ~ L C {N), these reunion probabilities, converge in a proper scaling limit, when 
N — >• 00, to a limiting form which can be expressed in terms of the Tracy- Widom distributions 
T\ and J^: this fact was certainly one of the main results obtained in Ref. [19]. Here, we have 
provided a detailed and self-contained derivation of these results. The main emphasis of the paper 
is on the study of the large deviations of these reunion probabilities both in the right and in the 
left tail, together with a careful analysis of the matching between the different regimes (left tail, 
central part and right tail) . While this matching in model A and model C is very similar to the one 
found for the distribution of the largest eigenvalue of GUE, one finds that the situation is much 
more involved in model B. In this case, there is instead a crossover (see Fig. [3]) between the central 
regime, for L — 2\/N ~ C(A -1 / 6 ) and the right tail, L > 2y/N, of the reunion probability, this 
crossover happening at a crossover length scale L cmss — 2VA cx A _1 / 6 (ln A) 2 / 3 , which seems to be 
a peculiar feature of these vicious walkers problems. 

Acknowledgements This research was partially supported by ANR grant 2011-BS04-013-01 WALKMAT 
and in part by the Indo-French Centre for the Promotion of Advanced Research under Project 4604-3. 



A Details about the large deviation regime 

In this appendix, we give some details concerning the calculation of Fjv(L) in the large deviation regime. 



A.l Derivation of the formula given in Eq. (|141D 

We first provide a derivation of the formula given in Eq. (|14ip starting from Q6ip and (|tJ7[) . Indeed, from 
the definition of Rk(a) one has 



hi(a) 



Ri(a)ho(a) , /12(a) = 7?2(a)/ii(a) = ho(a)Ri(a)R2(a) ■ ■ ■ 



(221) 



and more generally 



h k (a) = h (a)R 1 (a)R 2 (a) ■ ■ ■ R k {a) . 



(222) 
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Therefore, On (a) in Eq. (|96[1 can be rewritten as 

JV 

0(a, TV) = JV! J] h 2j -i(a) = AT! [Ao(a)«i(a)][fto(a)«i(a)ife(a)Ji8(o)] 
i=i 

x [/io(a)-Ri(a)iZ 2 (a)i?3(a)ii4(a)-Rs(a)] ■ • ■ 

= iV![/ l0 (a)i?i(a)] iv [7? 2 (a)i?3(Q)] iV - 1 [i?4(a)i? 5 (Q)] iv - 2 ■ • ■ 

JV-l 

= Nl[ho(a)Ri(a)] N J] [ J R2 fc (a) J R 2fc+1 (a)] iv - fe (223) 
fe=i 

which is the formula given in Eq. (|141[) . 

A. 2 Derivation of the formula given in Eq. (|142p 

In this appendix we give a detailed derivation, starting from the exact expression of Fn(L) in Eq. (|61II62[) . 
of the asymptotic estimate for lnFjv(L), valid in the limit L 3> y/2N 2> 1: 

In F N (L) = e" 5 [G 2N (£) + I 2N (£)] , (224) 

in terms of the variable £ = 7r 2 /a = 2L 2 and the functions G2iv(£) given in Eq. (|127[) and I2n(£) given by 

" 2fc + 1 



/ aw (0 = -2€$:%#- (225) 



This formula (|224[) was given in Ref. |31) [see their Eq. (28)] without any detail: that is the purpose of this 
appendix to fill this gap by providing a detailed derivation of it. 

To study the formula for Fn(L) starting from Eq. (|61l I62[) . we first write the product of the gamma 
functions in the denominator as 

II r(2 + j)r (| + /) = p k\ r Q + *) . (226) 

Note that the product of factorials can be written as 

N JV-l 

Y[k\ = 2 N ~ 1 ?, N - 2 ---N = N\ H k N ~ k . (227) 
fc=i fe=i 

Similarly, the product of F(k + 1/2) in Eq. (|226p can be written as 

n r ( k + 5) = f ^r 1 ) N n( fc + h N ~ k > ( 22§ ) 



where we have used 



2) V 2 y 1±v 2 



r(* + i/2) = (*-i)(*-|).-.ir(i) . (229) 



Finally, using Eqs (|227|l and (|228]) we write Eq. (|226jl as 

JV-l /o \ / /— \ N N-l 



Hr(2 + j)r(l+j)=N\(^-) U kN ' k U( k + b N - k ^ (230) 



3=0 
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where we have used T(l/2) = 0F. 

Using the ansatz for Rk(a) given in Eq. (|116|l where we consider Cfe(a)e _7r ^ a -C 1. with a = n 2 /2L 2 , 



we write, 



and for k > 1 



7T 1 



-R2fe(a)-R2fc+i(a) = 



ho(a)Ri(a) = J-— 1 + 2e~ 
a 2a 



fc A; + 1/2 r a 

l + -r c 2fc( a ) e 

a a L fe 



1 + 



1 + 2aci (a)e 



-c 2fc+ i(a)e 



(231) 



+ 0(e^ /a ) ,(232) 



fc + 1/2 

so that, using the formula in Eq. (|230|) one obtains that F?j(L) in Eq. (|61l I62[l can be written as 



F N {L) = l + 2e 



1 + 2qci (a)e 



(233) 



N-k N-l 

n( 



fc + 1/2 



c 2fe+ i(a)e 



2 \ N-k 



We now perform the expansion of this expression (|233[l . considering Ck(a)e~ n l a <C 1. It is more convenient 
to expand its logarithm ln_Fjv(L) which using the explicit expression for a (a) in Eq. (| 1 1T|) together with 
the expression of Ck(a) in terms of in Eq. (|120|l . can be written as a function of f = tv 2 /a as 



In JV(I) = -2N(2( - l)e~ e - lie-' £ ^iL-*o„(Q + JL-A G „ +l(S ) 



2fe + 1 



\fc=l fc=0 

where we have used Gi(£) = 1. It is then possible to write these two sums in Eq. (|234[) as 

lnFjv(L) = 2Ne-t - 2*r« ^ ^G,(€) + £ j . 

Finally, using the following identity satisfied by the polynomials Gk(£) [30] 

2iV-l 9N _ . 

2N -2£ ^GjiO = G 2N (0 , 

i=i 3 



(234) 



(235) 



(236) 



which can be shown, for instance, by using their explicit expression of Gfc(£) (| 12T|) . one obtains finally 

InFjv(L) = e" s (G 2JV (£) + h N (0) , (237) 

as given in the text in Eq. (|142|) . 



A. 3 An integral representation for I2n(£) in Eq. (11431) 

We start with the integral representation for Gfc(£) in Eq. (|126|) to express J 2 jv(£) in Eq. (|143|) as 

dt 



hniO = -2£<p 



-2fj 



Ci 



2ni 



S»+t( 1 + * 



2fc+l 



(238) 
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Performing an integration by part one obtains 



r „-2£< N ~ 1 / i \ 2fe 

^K^VEK) ■ (239) 



By performing the sum over k in (|239|) and drop ping an iV-independent constant term, one finally 
arrives at the expression given in the text in Eq. (|144[) . 
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